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Abstract. We investigate certain classes of integrable classical or quantum spin 
|' systems. The first class is characterized by the recursively defined property P saying 

' that the spin system consists of a single spin or can be decomposed into two uniformly 

coupled or disjoint subsystems with property P. For these systems the time evolution 
can be explicitcly calculated. The second class consists of spin systems where all 
non-zero coupling constants have the same strength (spin graphs) possessing N — 1 
independent, commuting constants of motion of Heisenberg type. These systems are 
' shown to have the above property P and can be characterized as spin graphs not 

. containing chains of length four. We completely enumerate and characterize all spin 

graphs up to N = 5 spins. Applications to the construction of symplectic numerical 
ly-j ■ integrators for non-integrable spin systems are briefly discussed. 

o ' 

Oh; 

^ (Jh 1. Introduction 

> 

Classical spin systems are examples of Hamiltonian mechanical systems. Hence the 
term "integrable" has a precise meaning in the context of the Liouville-Arnold theorem 
PP. It requires that there exist A" independent, commuting constants of motion, where 
A" denotes the number of spins and "commutation" is understood w. r. t. the Poisson 
bracket. For integrable systems one can find so-called action- angle variables I n , (p n such 
that 

r dH dH , r 

I n = -- — = 0, <p n = — =a> n = const. ,n = l,...,N . (1) 

dif n dl n 

Hence the equations of motion can be solved explicitely if the integrations involved in 
the definitions of the I n ,^Pn can be performed, see pQ. Since integrable systems are 
notoriously rare it is important to have as much examples as possible in order to test 
conjectures about larger classes of systems. In the last decades A^-soliton solutions of 
integrable infinite-dimensional Hamiltonian systems have found considerable interest, 
see e. g. These include solitary spin waves in infinite chains in the continuum limit 
[Sj . This motivates the search for other classes of integrable spin systems and for criteria 
of integrability. The task to independently characterize the class of all integrable spin 
systems XS seems extremely difficult. Most published work on integrable spin systems 
deals with special examples and numerical case studies jljjEj- Also in this article we will 
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not characterize IS itself, but certain subclasses of IS. 

For quantum systems the corresponding notion of "integrability" is less precise. 
Although some prominent examples of integrable classical systems have solvable 
quantum mechanical counterparts, including the harmonic oscillator, the Kepler 
problem and the two center Kepler problem, there is no comparable general theory 
of integrable quantum systems. For example, the general Heisenberg spin triangle with 
different coupling constants is an integrable classical system, but we do not know of 
any procedure to analytically calculating the eigenvectors and eigenvalues of the general 
quantum spin triangle for arbitrary individual spin quantum number s. However, for 
the subclasses of integrable spin systems to be considered below the eigenvalue problem 
of the corresponding quantum spin Hamiltonian can be analytically solved. 
In this article we investigate a subclass HIS C IS of the class of integrable spin 
systems, called Heisenberg integrable systems, or, shortly, H-integrable systems. They 
are defined by the extra condition that iV — 1 of the N constants of motion, as well as 
the Hamiltonian H itself, are of Heisenberg type, i. e. consist of linear combinations of 
scalar products of spin vectors: 

£M = 5><W*m ' ^ n = l,...,iV-l. (2) 

The remaining iV-th constant of motion is chosen as the 3— component of the total spin 
. We conjecture that HIS = IS if H is of Heisenberg type, but we have not proven 
this although we have some evidence from numerical studies of Ljapunov exponents 
for small spin systems, see also [Hj. We didn't obtain an independent characterization 
of HIS itself, but only for two subclasses HIQ C BS C HIS called "Heisenberg 
integrable spin graphs" and "^-partitioned spin systems". "Spin graphs" are systems 
with Heisenberg Hamiltonians 

H = ^2 J ^s» ■ $v , (3) 

fJ,<U 

satisfying 6 {0, 1}. Obviously, the coupling scheme of such a system can be repre- 
sented by an undirected graph, the N vertices of which correspond to the N spins and 
the edges (/i, v) to those pairs of spins where J^ u = 1. 

All spin graphs with N < 4 turn out to be H-integrable, with the exception of the 
4-chain. One main result of this article is that a spin graph is H-integrable iff it con- 
tains no 4-chain as a subsystem iff it is the union of two uniformly coupled or disjoint 
H-integrable subsystems ("uniform or disjoint union"). By recursively applying the uni- 
form union property we obtain a partition of the whole spin graph into smaller and 
smaller H-integrable subsystems with uniform or vanishing coupling. This sequence 
of partitions can be encoded in a binary "partition tree" B. Removing the condition 
Jp,v £ {0, 1} we arrive at the slightly more general notion of i3-partitioned spin systems 
for which the time evolution can be analytically calculated. The observation that the 
uniform union of two integrable systems is again integrable is certainly not new, see 
e. g. [TjjS] for a special case concerning quantum spin systems or for classical sys- 
tems. However, the use of partition trees in order to obtain the time evolution or the 
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eigenvectors in closed form seems to be novel. 

Our article is organized as follows. In section 2 we present the pertinent definitions 
and first results on classical integrable or H-integrable spin systems. Further results on 
subsystems and uniform unions are contained in section 3. Among these is theorem 1 
saying that any subsystem of an H-integrable spin system is again H-integrable. 
Section 4 contains our results on spin graphs. For example, theorem 2 states that each 
H-integrable spin graph is the uniform union of two H-integrable subsystems. Finally 
we will prove that a spin graph is H-integrable iff it does not contain any 4-chain, see 
theorem 3. As an application, we enumerate all connected spin graphs up to iV = 5 in 
the appendix. If they are H-integrable we indicate the uniform decomposition as well 
as the N commuting constants of motion; if they are not we display some 4-sub-chain. 
The next section 5 is devoted to the explicit form of the time evolution for £>-partitioned 
spin systems, see theorem 4. It turns out that their time evolution can be described 
by a suitable sequence of rotations about constant axes. This is closely related to the 
definition of action-angle variables satisfying (0), as we will show in section 5.2. 
In section 5.3 we will sketch how to calculate the eigenvectors and eigenvalues of the 
Hamiltonian in the quantum version of £>-partitioned spin systems. Section 6 contains 
a summary and an outlook. 

2. Definitions and first results 

Classical spin configurations are most conveniently represented by iV-tuples of unit 
vectors s = (s\, . . . , S/v), \s^\ 2 = 1 for /x = 1, . . . ,N. The compact manifold of all such 
configurations is the phase space of the spin system 

V = V N = {(s u ...,s N ) ||s M | 2 = 1 for fi= l,...,iV} . (4) 

The three components s* (i — 1, 2, 3) of the /i-th spin vector can be viewed as functions 
on V 

S ;-.V^R. (5) 

In order to formulate Hamilton's equation of motion we need the Poisson bracket 
between two arbitrary smooth functions 

f,g:V-^R. (6) 

The Poisson bracket has to satisfy a couple of general properties, namely bilinearity, 
antisymmetry, Jacobi identity and Leibniz' rule, see e. g. jTHj 10.1. Hence it suffices to 
define the Poisson bracket between functions of the form (J3J): 

3 

K,4} = ^2^jks^, (7) 

k=l 

where 5^ denotes the Kronecker symbol and the components of the totally 
antisymmetric Levi-Civita tensor. This definition turns V into a Poisson manifold. 
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We will sketch a more abstract way to endow Vn with a Hamiltonian structure: Let 
G be a Lie group and g* the dual of its Lie algebra, g* is endowed with a canonical 
Poisson bracket and, moreover, is a disjoint union of the orbits of the co-adjoint action 
of G upon g*. Every such orbit is a natural symplectic manifold, see [TU], chapter 14. 
The phase space Vn of a classical spin system results if G is taken as the iV-fold direct 
product the rotation group 50(3). In this case g* = (IR 3 )^ can be given the structure 
of a product of Euclidean spaces, unique up to an arbitrary positive factor, such that 
G operates isometrically on g* . Then Vn is the co-adjoint orbit consisting only of unit 
vector configurations. Hence the in (J7|) has its origin in the Lie bracket of the 
Lie algebra of 50(3), which is also the origin of the commutation relations of angular 
momenta in quantum mechanics. In the sequel we will, however, not make use of this 
abstract approach. 

Having defined the Poisson bracket, we can write down the differential equation 
corresponding to a given smooth function H : V — > R: 

±si = {s%H}, /i = l,...,iV, * = 1,2,3. (8) 

The r. h. s. of (JHJ) can be viewed as the vector field Xjj on V generated by the function 
H. If H is the Hamiltonian of the spin system, (jHJ) is Hamilton's equation of motion and 
Xh is the corresponding Hamiltonian vector field. It is complete since V is compact, 
see pi] Cor. 4.1.20. Generally, we define 

F t {H)s{0) = s(t), teR, (9) 

where s(t) = (si(t), . . . , sV(t)) is the solution of (jHJ) with initial value s(0). 

Ft{H) : V — ► V is called the flow of H and is defined for all t e R due to the 

completeness of the Hamiltonian vector field. 

Lemma 1 Let H,K : V — > R be smooth functions. Then the flows J~t{H) and J~t{K) 
commute iff {H, K} = 0. 

Proof: The if-part is a standard result, since the commutation of the flows is equivalent 
to = [X H ,X K ], see [HI 4.2.27, and [X H ,X K ] = -X {H)K}l see PH] 5.5.4. For the only- 
if-part we conclude = [Xh,Xk] = —X{h,k}, hence {H,K} = c = const.. For c / 
the differential equation 

j t K(s(t)) = {K,H} = -c (10) 
has unbounded solutions, which is impossible for compact V . Hence {H, K} = 0. □ 

For the rest of this section we consider a fixed Heisenberg Hamiltonian H : V — > R. 
It will be convenient to identify the spin system with its Hamiltonian. 
A constant of motion is a smooth function / : V — > R which commutes with the 
Hamiltonian: {/, H} = 0. if is said to be of Heisenberg type, or, short, a Heisenberg 
Hamiltonian, if it is of the form 

H(s) = ^ JfjySp ■ s v . (11) 
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The real numbers are called coupling constants. It will be convenient to set J ViX = 
for ji < v. Define the total spin vector 

N 

with components i = 1,2,3 and square S 2 = S ■ S. A Heisenberg Hamiltonian 
commutes with all components of the total spin and its square: 

= {H,S 2 } = {H,S {l) }, i = 1,2,3 . (13) 

Let A 7^ be any subset of {1, . . . , N} and a = \A\. Then Va denotes the phase 
space of the subsystem A, i. e. the manifold of all spin configurations of the form 
= (s M1 , . . . , s Mo ) such that /ii < /i 2 < . . . < fi a and |s Mi | 2 = 1 for all /!« G A. The 
Hamiltonian Ha of the subsystem A will be defined by 

Similarly, we define Sa = SugA together with its components and its square S\. 
Next we consider a decomposition of {1, ... , N} into two disjoint subsets, {1, . . . , N} = 
AJB such that A,B ^ 0. Let further N A = \A\ and N B = \B\. The Heisenberg 
Hamiltonian H is accordingly decomposed into three parts: 

+ 

= H a + H b + Hab ■ (16) 

Here and in the sequel we identify functions of Heisenberg type defined on the phase 
space of a subsystem Va with their unique extension to the total phase space V, defined 
by Jfj, u = for all /i, v with fi ^ A or v (jz A. 

If the coupling constants J^, fi G A,u G B occurring in Hab have all the same non-zero 
value, say, J^ u = cab G K, c^b 7^ 0, we will call the system if the uniform union of 
the subsystems Ha and H B . If the analogous condition holds with cab = we call the 
system H the disjoint union of the subsystems H A and H B . A Heisenberg system is 
called connected if it is not the disjoint union of two subsystems. A Heisenberg system 
where all coupling constants are non-zero and have the same value will be called a 
pantahedron. 



Definition 1 A Heisenberg spin system H is called Heisenberg integrable, or, short, 
H-integrable, if there exist N — 1 independent constants of motion E^ of Heisenberg 
type which commute pairwise: 

{E {n \E (m) } = for alln,m = 1,...,N-1 . (17) 

The above condition of independence means that there exists some s G V such 
that the set of covectors {dE^(s), . . . ,dE^ N ^(s)} is linearly independent. It follows 
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that this condition is also satisfied in some neighborhood of s G V . But it cannot 
hold globally: If you take some F in the linear span of the E^ such that s G V is a 
critical point of F, i. e. dF(s) = 0, then it is obviously violated. Later we will derive 
a simple criterion for the independence of a number of Heisenberg constants of motion, 
see proposition [T] 

If a connected spin system is H-integrable, it can be easily shown that 
{E^\ . . . , E^' 1 ^, E^ = S^} will be a set of N independent, commuting constants of 
motion. Hence any H-integrable system is also integrable in the sense of the Liouville- 
Arnold theorem. We conjecture that the converse is also true. 

For an integrable spin system the N constants of motion E^ are not uniquely 
determined. First, one can consider linear transformations 

N 

F^ = J2^nmE im \ n = l,...,N, (18) 

m=l 

where the A nm are the entries of an invertible matrix. These transformations leave in- 
variant the space £ of functions spanned by the E^ n — 1, . . . , N. But also the space £ 
need not be uniquely determined by the Hamiltonian H: Consider the disjoint union H 
of two H-integrable subsystems Ha and Hb- Then one could either consider the union 
of the two sets of independent, commuting constants of motion for the subsystems, in- 
eluding S A and S B , or, alternatively, the union of the Heisenberg constants of motion 
of the subsystems together with S 2 — S\ — S B and S^ 3 \ Since the second choice of the 
E^ is always possible, our definition ^ entails that the disjoint union of H-integrable 
spin systems will be again H-integrable, see proposition El 

Lemma 2 For any spin system with N spins there exist at most N independent, 
commuting constants of motion. 

Proof: Consider a set of M independent, commuting constants of motion E^ n \ 
not necessarily of Heisenberg type, and the condition of independence holding in some 
neighborhood of s G V. Then the equations E^ n \s) = E^(s ), n = 1, . . . ,M define 
a local 2N — M-dimensional submanifold S of V and the vectorfields X E („) are linearly 
independent and tangent to the submanifold S. It follows that M < 2N — M, i. e. there 
exist at most iV independent, commuting constants of motion. □ 

A spin graph is a Heisenberg spin system where all non-zero coupling constants 
have the same strength J ^ 0. Without loss of generality we may assume J^ u G {0, 1}. 
As explained in the introduction, the system can be represented as an undirected graph 
with vertices. Of course, the above definition of a connected Heisenberg spin system 
coincides with the graph-theoretic notion of connectedness if the system is a spin graph. 
The notion of a subsystem used in this article is, in the case of spin graphs, equivalent 
to what graph theorists call a "vertex-induced subgraph", see [12]. It means that the 
subsystem is obtained by removing a number of vertices along with any edges which 
contain a removed vertex. The following fact is well-known, see [T2"] . theorem 1.6: 
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Lemma 3 In a connected spin graph with N > 2 vertices one can remove two suitable 
vertices such that the remaining subsystem is still connected. 

In order to evaluate the Poisson bracket between functions of Heisenberg type and to 
argue with the resulting equations we need the following lemmas which are easily proven: 

Lemma 4 (i) {s M • s u , s x ■ s u } = ■ (s x x s„) = det(s /i , s x , s u ) 
Let k,£ £R 3 be constant vectors, then 

(ii) {s M ■ s u ,k- s u ] = det(s At , k, s„) 

(iii) {k ■ s u , £ ■ s u } = det(k, £, s v ) 

Lemma 5 If one of the following equations 

N 



E ^ • s» = 0, (19) 

N 

C^s^-s u = 0, (20) 

fl < V 

N 

■ {s„ x s u ) = 0, (21) 

ft < v 
N 

C^ x det{s fl ,s u ,s x ) = (22) 

fi < v < A 

holds for all (si, . . . , sn) G V then all coefficients of the corresponding equation must 
vanish. 

Proof: By induction over iV using the replacement sjy+i )—* — Sjv+i- 

□ 

Now we can formulate a criterion for the commutation of two functions of 
Heisenberg type: 

Lemma 6 Let H be a Heisenberg system and E a function of Heisenberg type. E 
commutes with H iff 

EfiviJftx — J vX ) + E flX (J uX — J^ v ) + E uX (J fll/ — J mA ) = (23) 

for all fj, < v < A < N. 

Proof: The lemma is proven in a straight forward manner by using lemma Efi), 
cyclic permutations of triple products and (|22p. □ 



Next we will show that the independence of M functions of Heisenberg type is 
equivalent to the linear independence of the corresponding symmetric matrices. 
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Proposition 1 Let E^ 1 ' : V — > WL be M functions of Heisenberg type, i. e. 

(b) = •* v ,n = l,...,M. (24) 

and denote by E^ n ^ the symmetric N x N -matrix with entries E^" = E$ and = 
for fj,, v — 1, . . . , N — 1. T/ien i/ie following conditions are equivalent: 

(i) There exists an s E V such that the set of covectors {dE^ n \s)\n = 1,...,M} zs 
linearly independent. 

(ii) T/ie set o/ matrices {E^|n = 1, . . . , M} is linearly independent. 

Proof: We will prove the equivalence of the negations of (i) and (ii). 

(i) =>- (ii). Assume that {E^ \n — 1, . . . , M} is linearly dependent. That is, there exists a 
non-vanishing real vector (Ai, . . . , Am) such that J2 n ^nM = 0. It follows that E(s) = 
£„A n £( n )(s) = for all s E P and hence = d (£ n KE in) (s)) = J2 n X n dE^(s). 
Hence the set of covectors {dE^ n \s)\n = 1, . . . , M} is linearly dependent for all s E V. 

(ii) =>■ (i). It is possible to invert the sequence of arguments of the first part of the 

? 

proof, except at the step dE = =>E = 0. Here we can only conclude E = c = const, 
and obtain the apparently weaker condition 

c = E(s) = An^ (n) (s) = X nE$s, ■ s v (25) 

for all s E V. Replacing s v by —s v for fixed v yields ^^E^s^ • s v = 0. By summing 
over v we obtain 5^ M<1/ E^s^ ■ s v — and, by (f2T)|) , £^ = for all fi < u. Hence 
]T n A n E( n ) = and the set {E^\n = 1, . . . , M} is linearly dependent. □ 

We note that (J2SJ) can be viewed as a system of (^) linear equations for the ( 2 ) 
unknowns E^ u . An H- integrable system admits at least N — 1 linearly independent 
solutions, according to proposition d In general, it will admit more solutions, but only 
iV — 1 of these will lead to commuting constants of motion. Hence, in the H-integrable 
case, the matrix M of the system of linear equations (|23j) has the rank r < ( 2 ) — (iV — 1). 
For N = A the rank condition r = („) — (N — 1) = 3 is even sufficient since it implies 
the existence of a constant of motion E which is not of the form XH + fiS 2 . After some 
algebra we obtain: 



Proposition 2 A Heisenberg system with N = 4 spins is H-integrable iff 

. (26) 



Jl3 ~ J23 J23 — J\2 
^14 — <^24 J24 — J\2 

J14 — J34 J34 — J13 



This criterion can be used to independently check the results on spin graphs with N = 4, 
see appendix A. 
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3. Subsystems and uniform union 

In this section we will collect some general results on H-integrable systems in connection 
with subsystems and uniform unions. 

Lemma 7 Let E be a Heisenberg constant of motion of a Heisenberg system H and 
0^i4c{l,..., N}. Then the restriction Ea is a Heisenberg constant of motion of the 
subsystem Ha- 

Proof: The lemma follows from lemma EH since the restricted functions Ea and Ha 
commute iff the equations (j2Sj) hold for \i < v < X with /z, v, A e A. □ 

Proposition 3 If a Heisenberg system H is the uniform or disjoint union of two H- 
integrable subsystems Ha and Hb, then H itself is H-integrable. 

Corollary 1 Each pantahedron is H-integrable. 

Proof: Let Ea be one of the Na — 1 independent, commuting Heisenberg constants of 
motion of Ha- In particular, Ea commutes with Ha and also with Hb since B n A = 0. 
Since Ea is a function of Heisenberg type it commutes with S 2 , S\ and S B , hence also 
with Hab since 

Hab = \cab (S 2 -S\- Si) - (27) 

It follows that Ea commutes with H . The same holds for a corresponding constant of 
motion Eb of the second subsystem. Hence the Na — 1 functions Ea together with the 
Nb — 1 functions E B and S 2 — S\ — S B form a set of N — 1 independent, commuting 
constants of motion of Heisenberg type. This means that H is H-integrable. □ 

The converse of proposition El is not true: The general spin triangle is a Heisenberg 
spin system with N = 3 and three different coupling constants. It has H and S 2 as 
independent, commuting constants of motion and is hence H-integrable, but it is not 
the uniform or disjoint union of two H-integrable subsystems. 
Next we will show that H-integrability is heritable to subsystems. 

Theorem 1 Any subsystem Ha of an H-integrable system is itself H-integrable. 

Proof: Consider N — 1 independent, commuting constants of motion Ex, ... , E^-x of 
the form 

E i = Y,E i ,^s ll -s u . (28) 
These constants of motions span a linear space T. 

According to the assumptions of the theorem {1, . . . , N} is the disjoint union of two 
nonempty subsets A and B such that N = Na + Nb with Na = \A\ and Nb = \B\. 
We arrange the coefficients Ei^ v in the form of a matrix E with ( 2 ) rows and N — 1 
columns. The first rows contain the coefficients with fi < v and \i, v 6 A; the next 
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( 7 2 B ) rows contain the coefficients with \x < v and //, v G B and finally the remaining 
NaN b rows those with fj, E A, u E B. In this way the matrix is divided into three 
blocks, see the following figure. 

Ei E2 ... En-i 

]1 ( " A) 

• • • r 

N A N B 

N- 1 

Next this matrix will be transformed into a lower triangular form by elementary 
Gauss transformations. We allow arbitrary permutations of columns and arbitrary per- 
mutations of rows within the three blocks, see the following figure. 

Fi F2 ... Fn-i 

(*) 

N 4 N„ 







d B 


d A B 



N- 1 

The resulting matrix F begins with (1a linearly independent columns spanning a 
linear space J- a- <^a is the maximal number of independent constants of motions of the 
subsystem A obtained as restrictions to A of functions from T . The next ds columns 
of F span the linear space Tb- ds is the maximal number of independent constants of 
motions of the subsystem B obtained as restrictions to B of functions from T which 
vanish on A. The remaining dAB columns span the linear space Tab of functions from 
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T vanishing on A and B. Since elementary Gauss transformations do not change the 
rank of the matrix, we have 

d A + d B + d AB = N -1 . (29) 

For sake of simplicity we identify the columns of F with the corresponding constants 
of motion. By lemma the restrictions to A of Fi, . . . , F dA will be constants of motion 
of H A . According to lemmaElthere are at most N A — 1 independent constants of motion. 
The analogous argument holds for F dA+1 , . . . , F dA+dg and the subsystem H B . Hence 

d A < N A - 1 and d B < N B - 1 . (30) 

It follows that N - 1 = d A + d B + d A B <N A -1 + N B -1 + d AB = N - 2 + d AB , 
hence 

d AB > 1 • (31) 

Next we want to show that d AB < 1. In this case we are done: d A < N A — 1 would 
imply N — 1 = d A + d B + d AB < N A — 1 + N B — l + l = iV — 1 which is a contradiction. 
Hence d A = N A — 1 and the subsystem H A would be H- integrable. 



Proving d AB < 1 is equivalent to show that T AB is at most one- dimensional, i.e. 
that the ratios of the coefficients F^\/F UK with G A and A, k G B are uniquely 
determined. Hence consider some F G {F dA+dB+ i, . . . , F N _i}. Its restrictions are 
F A = F B = 0. For v G A and A G B lemma |B1 yields 

F^xiJuX — Jftu) + F U \(J^ U — J^x) = 0. (32) 
Then the following ratio of coefficients 



Ffj, \ J jj,v J^iX 



(33) 



F v x J V X J fjLU 

is uniquely determined, except in the case where the nominator J^ v — J^x and the 
denominator J v x~ J^v vanish simultaneously: 

J^u = J^x = Jux ■ (34) 

Define, for fixed /i G A and A G B, M(yU, A) to be the set of v G A satisfying (J34)) . If 
M(fi, A) = A then H A is a pantahedron and hence H-integrable. If M(/i, A) 7^ A then 
there exists a k G A such that k ^ M(/x, A) and the ratios F^x / F kX and F uX / F kX are 
uniquely determined. Hence also the ratio 

Ffj,x F^x F K x (35) 
F v \ F K x F v x 

will be uniquely determined. By analogous reasoning, also the ratio F v x / F VK with 
v G A and A, n G B, and hence F^x/F VK with \i,v G A and A, n G B will be uniquely 
determined. This completes the proof. □ 



The preceding proof also shows: 
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Corollary 2 Let H be an H-integrable spin system with two subsystems Ha and Hb 
such that {1, . . . , N} = AOB. According to theorem^ Ha and Hb are also H-integrable. 
Then the N — l independent, commuting constants of motion Fx, ... , can be chosen 
such that 

(i) Fi,...,Fn a -i are independent, commuting constants of motion of Ha and vanish 
on B, 

(ii) Fn a , ■ ■ ■ ,Fn~2 o,re independent, commuting constants of motion of Hb and vanish 
on A, 

(Hi) Fjv-i vanishes on A and on B. 

Proof: It remains to show that the Fx, . . . , F^ A _x of (i) vanish on B. This can be done 
by further Gauss transformations which add suitable multiples of columns of (ii) to the 
columns of (i). □ 



4. Spin graphs 

As explained in section 2, spin graphs are Heisenberg spin systems such that the coupling 
constants satisfy J^ v G {0,1}. For these systems, H-integrability can be completely 
analyzed. According to proposition El the uniform or disjoint union of H-integrable 
systems is again H-integrable, but not all H-integrable systems are obtained in this way. 
However, all H-integrable spin graphs are the uniform or disjoint unions of H-integrable 
subsystems, as we will show. This means that there is a construction procedure by 
which one can compose all H-integrable spin graphs from small constituents. Starting 
with two single spins, which are trivially H-integrable, we can either form a disjoint 
union or a spin dimer. The uniform union of a dimer and a single spin yields a uniform 
triangle; the uniform union of a single spin with a pair of disjoint spins is a 3-chain, etc.. 
Remarkably, the 4-chain cannot be obtained in this way and is hence not H-integrable. 
Our first result is: 

Lemma 8 Each connected H-integrable spin graph with N > 1 vertices is the uniform 
union of two subsystems. 

Proof: The proof will be performed by induction over N. For N = 2 the theorem 
holds since the dimer is the uniform union of two single spins. 

Next we assume the theorem to hold for all spin graphs with iV or less vertices and 
consider a connected spin graph with N + 1 vertices and Hamiltonian H. According to 
lemma El we may assume that the subsystem H N with vertices {1, . . . , N} is connected 
and, by theorem Q H-integrable. We denote by H^+x the single spin system with vertex 
N + 1. By the induction hypothesis and since N > 1, H^ is the uniform union of two 
H-integrable subsystems Ha and Hb, where {1, . . . , N} = AijB and A, B ^ 0. 
Ha is further decomposed into subsystems H A and H\ with vertex sets A and Ax, 
respectively. A consists of those spins in A which do not couple to iV + 1; Ax consists 
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of the remaining spins in A which thus uniformly couple to iV + 1. The analogous de- 
composition is performed w. r. t. Hb- 

Since A, B ^ we must have Aq ^ or A\ ^ 0, and, similarly, Bq ^ or B\ ^ 0. In the 
case A = B = the proof is done since this means that H N couples uniformly with 
H^+i- The case Ax — B\ — can be excluded since it implies that H is disconnected. 
Hence it suffices to consider the case Ai ^ and Bq ^ in what follows. The other 
remaining case A ^ and Bi ^ can be treated completely analogously. The situation 
is illustrated in the following figure. 




If Aq = 0, the total system if is a uniform union of the two subsystems with 
vertex sets A and {A^ + 1}U5 and the proof is done. Hence we may assume A ^ 0. 
If the coupling between the subsystems H\ and H\ is uniform we may rearrange the 
decomposition by setting A' = and B' = Aq U Bq, leaving Ai and B\ unchanged. 
But this case has already be considered above. Hence we may assume that the coupling 
between H\ and H\ is non-uniform. 

By corollary El (iii) and since is H- integrable, the total system H possesses a non-zero 
constant of motion of the form 

E = ^ e ^n+i ■ s N+1 . (36) 

li<N+l 

Lemma El implies 

E^N+l{J^\ — J\N+l) + E\M+l(Jfj,N+l — JfMX) = (37) 

for all fi < A < N + 1 since E^\ = 0. We will show that the case A , At, Bq ^ is in 
contradiction to the above-mentioned fact that E is non-zero. To this end we consider 
(|3~7jl in the following four cases: 



jiQ G Bq and Ai G Ai (J, 



Mo Ai 



Ai N+l 



1 and J, 



Mo N+l 



0) 



E\ t N+l 



. 



(38) 



fix G B x and A G Aq (J Mi a = J^N+i = 1 and J Xq n+i = 0) 



=>• n + i — 



(39) 
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G B Q and A G A (J Mo jv+i = </a tv+i = and J t 



Mo -^o 



A G A , Ai G Ai and J AoAl = (J Ao jv+i = and J Ai tv+i = 1) 



(40) 



^E XoN+1 = 0. (41) 

Since the coupling between H\ and H\ is non-uniform, there exist A G A and 
A' : G Ai such that J\' Q x' 1 — 0. For this A the coefficient Ey n+i vanishes by (JHJ). 
Equation (f4T)|) then yields £? M0 at + i = for all /io G -Bo and, further, i£ Ao at + i = for all 
A E A . By the equations (f3H)l and the remaining coefficients of vanish, which 
leads to a contradiction. □ 



Lemma IHl theorem ^ and proposition El together yield: 

Theorem 2 Each H-integrable spin graph is the uniform or disjoint union of two H- 
integrable subgraphs. 

It follows from theorem El that all spin graphs with N < 3 are H-integrable, but 
that the chain with N = 4 is not H-integrable since it is not the uniform union of 
smaller systems. By virtue of theorem every spin graph containing a 4-chain will 
not be H-integrable. The converse is also true and yields the following graph-theoretic 
characterization of H-integrable spin graphs. 

Theorem 3 A spin graph is not H-integrable iff it contains a chain of length four. 

Proof: It remains to show the only-if-part. This will be done by induction over N. 
For N = 4 the theorem is proven by a complete classification of all connected spin 
graphs, see the appendix. Next we assume that the theorem holds for all spin graphs 
with iV or less spins and consider a spin graph H with N + l vertices which is not 
H-integrable. 

If H is the union of two disjoint subsystems, at least one of them is not H-integrable by 
proposition El and hence contains a 4-chain by the induction assumption. Thus we may 
assume that H is connected. By virtue of lemmaElwe can remove a suitable vertex with 
number, say, N + 1, such that the remaining subsystem is still connected. Further 
we may assume that is H-integrable, since otherwise it would contain a 4-chain by 
the induction assumption. The coupling between and H^ + i is neither uniform nor 
zero, since then H would be integrable by proposition El or disconnected. Hence we may 
decompose H N into a maximal subsystem H% which is not coupled with H N+ i and the 
remainder H\j which is uniformly coupled with Hjf+i- Both subsystems are non-empty 
and H-integrable by theorem ^ 
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H 




N+1 



if TV is connected and H-integrable and hence, by theorem |21 the uniform union of 
two non-empty subsystems Ha and Hb- Both subsystems are further decomposed into 
H A , H\ and H B , H B according to their coupling with H N+ i, similarly as H N above. 
Let A , Ax, B , B\ be the corresponding subsets of {1, . . . , N}. A U B ^ and 
A\ U B\ 7^ 0, see above. 




We consider the case A = 0. This means that Hn+i as well as H b is uniformly coupled 
to Ha- The restriction of H to {iV + 1} US cannot be H-integrable, since then H would 
be H-integrable by proposition El Hence {N + 1}UB contains a 4-chain by the induction 
assumption. Analogously we can argue in the case B = 0. 
Hence we may assume A ^ and B ^ 0. 

Assume that Hn is a pantahedron. Then H would be the uniform union of H]^ and the 
disjoint union of H N+ i and H% and hence H-integrable by proposition El contrary to 
previous assumptions. 




Thus H is not a pantahedron and possesses at least one uncoupled pair of spins, say 
(/!, v) such that J^ u = 0. We have /i G A , v G A\ or /i G B , v G Bi, since A and 
i? are uniformly coupled. Because of A ^ and -Bo ^ we can choose any A G A or 
A G Bq an d obtain a 4-chain (iV + 1, u, A, /x). 
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• > 

Ha" 










H B 



□ 



5. Time evolution of ^-partitioned systems 

5.1. Explicit form 

In this section we consider spin systems which are the uniform or disjoint union of 
subsystems A and B in such a way that these subsystems enjoy the same property, and 
so on, until one or both subsystems consist of single spins. Thus we obtain a nested 
system of partitions which can be encoded in a binary partition tree B. Examples are 
H-integrable spin graphs which are special £>-partitioned systems by virtue of theorem 

121 

The time evolution of such systems can be exactly described by means of a recursive 
procedure, see In this section we will, additionally, provide an explicit formula for the 
general solution of the equations of motion which was assumed to be "too cumbersome" 
in using the notion of a "partition tree". 



Definition 2 A partition tree B over a finite set {1,...,N} is a set of subsets of 
{1, . . . , N} satisfying 

(i) ®£B and {l,...,N}eB, 

(ii) for all M,M' e B either M n M' = or M C W or M' C M , 

(in) for all M e B with \M\ > 1 there exist Mi, M 2 E B such that M = M 1 UM 2 . 

It follows from|2fii) that the subsets M l5 M 2 satisfying M = Mi(jM 2 in definition Hfiii) 
are unique, up to their order. Mi,M 2 are hence defined for all M e B with \M\ > 1. 
Mi and M 2 denote the two uniquely determined "branches" starting from M. It follows 
that B is a binary tree with the root {1, . . . ,N} and singletons {/i} as leaves. More 
general partitions into k disjoint subsets can be reduced to subsequent binary partitions 
and hence need not be considered. For all M e B there is a unique path 

V M (B) = {M f G B | M c M'} (42) 

joining M with the root of B. It is linearly ordered since M C M' and M C M" imply 
M' C M" or M" C M' by definition El (ii) . Especially, every element /i e {1, . . .,N} 
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belongs to a unique, linearly ordered construction path 

V„(B) = {M G B | y. G M} . (43) 

Although there are many different partition trees over a fixed finite set, all have 
the same size which can be easily determined by induction over N: 

Proposition 4 \B\ = 2N — 1 . 

For {1, . . . , N} we will denote by M the "successor" of M, that is, the smallest 
element of Vm(B) except M itself. To simplify later definitions we set {1, . . . , iV} = 
and denote by B the class of all successors, i. e. 

B = {M\MeB}. (44) 

It follows that \B\ = N. Later we will use B as an index set for iV action variables. 
For n 7^ v G {1, . . . , N} let G B denote the smallest set of B such that /i, v G M^, 
i. e. M^ u G B is the set where both construction paths of \i and v meet the first time. 
Consider real functions J defined on a partition tree 

J : B — ► R . (45) 

Then 

H = Y,J(M,u)s ll -s v (46) 

defines a Heisenberg Hamiltonian. The corresponding spin system will be called 
a ^-partitioned system or sometimes, more precisely, a (B,J)— system. The iV- 
pantahedron is a (B, J)-system where B is arbitrary and J is a constant function. As 
mentioned before, all H-integrable spin graphs are (B, J)-systems. For example, the 
spin square is obtained by the partition tree 

B = {{1, 2, 3, 4}, {1, 3}, {2, 4}, {1}, {2}, {3}, {4}} (47) 

and the function J with J({1, 2, 3, 4}) = 1 and J(M) = else. 

Recall that S M denotes the total spin vector of the subsystem M C {1,...,N} 
with length Sm and that Ft{H) denotes the flow map describing time evolution of spins 
according to a Hamiltonian H . It follows by lemma |H1 that the squares of the total spins 
corresponding to a partition tree commute: 

{S M , S 2 M ,} = for all M, M' <E B . (48) 

We consider rotations in 3-dimensional spin space: 

Definition 3 D(cD, t) will denote the rotation matrix with axis Q and angle \uj\ t. D(0, t) 
equals the identity matrix I. 

The proof of the following proposition is easy and will be omitted: 
Proposition 5 Let M, M' C {1, . . . , JV}. Then 

(i) Ft^S 2 M )=B(SMit), 

(ii) B0 M ; t)S M = Sm if M' C M or M' fl M = . 
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A short calculation shows that the Hamiltonian (|46|) of a (£>, J)-system can also be 
written as 

H = \ E J ( M ) s * - s k - $k) = \ E ( J ( M ) - > («) 

if we set J({fj}) = J(0) = for aU // G {1, . . . , JV}. It follows that 

T t {E) = J] F t (~(J(M)-J(M))S M ) = J] 0(S M , (J(M)-J(M))t) ,(50) 
AfeB MeB 
where the product is understood as the composition of flow maps or rotations and the 
order of the composition does not matter according to (J48j) and lemma ^ 
In order to calculate !F t {H)s^ we need only consider those factors in the product (|5U|) 
where fi G M, i. e. M G V^B). This follows by proposition E^ii). Moreover, if we 
choose a decreasing sequence of sets from left to right in the product (|5U|h we can write 
0(S M ,(J(M) - J(M))t) = ©(SW(0),(J(M) - J(M))t). Hence we have proven the 
following: 

Theorem 4 Lei H be a (B, J) -system. Then its time evolution can be written in the 
form 

s ~^) ^ILe^s) D (^(°)> ( J ( M ) - ^W)t)^(O) , (51) 

where the arrow above the product symbol denotes a product according to a decreasing 
sequence of sets M G V^(B) from left to right. 

We note without proof that the time evolution in the presence of a Zeeman term in 
the Hamiltonian of the form H + B ■ S , where B is the dimensionless magnetic field, is 
obtained by multiplying (Jo"Tj) from the left with 3(B,t). 

5.2. Action-angle variables 

Although theorem 0] gives a final answer to the problem of time evolution of a (£>, J)- 
system, it will be yet instructive to relate this result to the general form of the time 
evolution in terms of action-angle variables. We will give the result for the general case 
including a Zeeman term B ■ S = Be - S in the Hamiltonian but without detailed proofs. 
We write Sm = Sm^m for MgB. For M — {1, . . . , iV} = we set e = e and S = e- S. 
Tt{\S 2 M ) = B(§ M ,t)_Md {f,\S\ { } = {f,S M }S M imply T t {S M ) = B(e M ,t). Thus the 
functions Sm, M G B, generate rotations about the axes cm with unit angular velocity. 
Moreover, due to proposition El they represent iV commuting constants of motion and 
hence are good candidates for action variables. 

It remains to define suitable angles which change with unit angular velocity only 
under rotations about the axes cm- To this end consider the three unit vectors 
c*m, &Mi an d c-Mx- We may assume ejj ^ ±e M , and cm 7^ ±4fi since the action- angle 
variables need only be defined on an open dense subset of the phase space. Defining 
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and 

e v = e M x e e (53) 

we obtain an orthonormal frame (eQ,e v , e^) depending on the chosen point in phase 
space. We define polar and azimuthal angles 9m, fM by expanding cm x into this basis: 

e Ml = sin 9 M cos (fM^e + sin 8 M sin(p M e v + cos 8 M e M ■ (54) 

Then the following result holds: 

Proposition 6 The functions Sm,^Pm, M E B, define action-angle variables of a 
(B, J) -system. Especially, they satisfy 

Wm>, S m } = 5 M , M for all M, M' e~B . (55) 
5.3. The quantum case 

A given (B, J)-system and a spin quantum number s 6 {|, 1, §, . . .} uniquely specifies 
the corresponding spin system. All functions on phase space considered above have 
unique representations as Hermitean operators acting upon a (2s + 1 ^-dimensional 
Hilbert space H. To avoid complications we postulate a correspondence between 
functions on phase space and operators only for those functions which are sums of 
monomials of s 1 ^ with Poisson-commuting factors. Usually, this correspondence will be 
denoted by a sub-tilde, e. g. Sm, ,MeB. The commutator of operators corresponds 
to the Poisson bracket of the corresponding functions, setting H = 1: 

[f,g] = i{f L g}. (56) 

Hence for ^-partitioned systems the N operators S 2 M ,M e B, commute too. Together 

with = e ■ S, where e is the unit vector into the direction of the magnetic field, 
these operators constitute a complete system, that is, their common eigenvectors are 
uniquely determined by the corresponding eigenvalues 

S 2 M ® = S M (S M + 1)$ , (57) 

g(3) $ = £(3) $ ; ( 5 g) 

and span 7i. Hence we may write 

*=I(SmW- (59) 
Since the Hamiltonian is a function of the S 2 M , S (3) , 

H = \ £ (J(M) - J(M))S 2 m + BS® , (60) 
Me B 

the vectors (JST?|) are also eigenvectors of H and the corresponding eigenvalues can be 
read off from (joTIj) . 

The eigensystem {$} can be obtained by an multiple tensor product representation of 
SU (2) starting from the (2s + l)-dimensional irreducible representations corresponding 
to single spins and following the sequential partition of {!,... ,N} given by B. The 
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product representation of two irreducible representations corresponding to quantum 
numbers S±, S2 is spanned by vectors of the form 

\S,S {3) )= CG(S u S 2 ,S[ 3 \si 3) ;S,S®)\S u s[ 3) )® \S 2 ,S?) , (61) 

s[ 3) , s^ 3) 

where CG(. . .) denotes the Clebsch-Gordon coefficients, see, for example, jTHj 27.9.1. 
In order to write the product representation in compact form we introduce the notations 

B = {M e B I M ^ {1, . . . , N}} , (62) 
B x = {M G B I \M\ > 1} . (63) 

Then <3> can be written as 

<&=I(Sm)m g B> (64) 

— I Y\. CG(Smi, Sm 2 i ^Mxi ^m 2 'i ^Mi £>m) I l^i \ ■ ■ ■ > &n ) -(65) 

rtf (3)s „ \Me Bi / 



Here \S[ , . . . , 5^) denotes the product state which is a common eigenvector of all 
with eigenvalues S^, \i — 1, . . . , N. 

This concludes the diagonalization of (£>, J)-systems in the quantum case. 



6. Summary and outlook 

We have completely characterized H-integrable spin graphs by the graph theoretical 
property not to contain 4-chains. For these spin graphs as well as for the larger class of 
^-partitioned systems the time evolution can be explicitely calculated as a product of 
certain rotations about constant axes. But obviously this property is very rare. It means 
that four connected spins must not form a chain but rather close to build a triangle or 
a square. Hence this special type of integrability will either be satisfied for the majority 
of small spin graphs consisting of, say, four or five spins, see the appendix. Or it may 
be satisfied for larger systems which are close to the pantahedron type, where each spin 
is uniformly coupled to each other one. As a rule, spin lattices do not fall into this 
class. Possible direct applications of our theory will thus be confined to small clusters 
or magnetic molecules. 

There is, however, another possible application of H-integrable spin graphs to the 
construction of numerical integrators for non-integrable spin systems. Numerical 
integrators can be viewed as computable transformations of phase space which 
approximate the exact time evolution. It appears that those transformations which 
respect the symplectic structure of phase space or, equivalently, its Poisson bracket, lead 
to numerical integrators with favorable properties, see ^3] . They are called "symplectic 
integrators". The splitting of a non-integrable Heisenberg Hamiltonian into integrable 
parts gives rise to symplectic integrators by means of Suzuki- Trotter decompositions, 
which exist of various orders, see [To] . Such a splitting is always possible, just consider 
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the decomposition into dimer Hamiltonians. But we expect that the numerical integrator 
is the better the larger the parts of the splitting are. In this situation the present results 
on H-integrable spin graphs including the explicit form of the time evolution can be 
utilized to develop efficient symplectic integrators for non-integrable spin systems. We 
have obtained first results in this direction which will be published elsewhere. 
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Appendix A. Connected spin graphs with N < 5 

Table lA"Tl contains all connected spin graphs up to iV = 5 spins. We have given a maximal number of 
independent, commuting constants of motion, using the abbreviation s?- = Sj • Sj. For H- integrable 
spin graphs a decomposition into uniformly coupled subsystem is indicated. If the system is not 
H-integrable the decomposition shows a sub-chain of length 4. 



N 


System 


Constants of motion 


H-integrable 


Decomposition 


1 


• 




yes 




2 




H, 


yes 


•'o 2 


3 






yes 


• o 2 m 3 




l 

• 


u o(3) 2 

H, SW, sf 3 


yes 


1 

f 

Wo 2 


A 


4 1 

ft 


tt o(3) „2 2 
) D ) s 13> S 24 


yes 


4 1 

• 

b m 2 




u 




yes 


• 






H, S^, 4, ^ 


yes 










yes 


• 

• 4 o 2 • 




f 






• 




A 


if, ^ 3 ), S| 4 , S?3 + S? 4 


yes 


o 2 








no 






a 1 ^ 2 ^ 3 w 
• • • • 


rr c(3) 2 _i_ 2 _i_ 2 

fi, & , s 13 + s 14 + s 24 






5 


b 


If, S 1 ^- 1 , sf 3 , S24 + sf 5 + S25 + ^35 


no 


b 




b 


H, S^ 3 \ sf 3 , S24 + sf 5 + S25 + S35 


no 


b 














b 


H, S^ 3 \ sf 3 , s\ 2 + s| 3 ,sf 5 + S25 + S35 


yes 


b 
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1 2 3 4 5 

• A • 

X 



A 
V; 

X 
X 

A 



■& 



TJ .9(3) s 2 s 2 4- s 2 s 2 4- s 2 4- s 2 



/J, 5( 3 ), S 2 3 + S? 4 + S 2 5 + s| 4 + S 2 5 + S, 



35 



/J, 5^, s 2 3 + s 2 u + s 2 15 + s| 5 + s 



35 



/J, 5( 3 ), S 2 3 + S 2 4 + S 2 5 + S 2 4 + S 2 5 



if, S( 3 ), S 2 3 + S 2 4 + 



S 15 + S 24 + S 25 



if, S-( 3 ), S 2 3 + S 2 4 + ^4 + si 5 + S 



35 



H, S^ 3 \ s 2 2 , s§ 4 , s 2 3 4- s 2 4 + s 23 4- s 
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if, S**- 3 ^, s 2 2 , s| 4 , s 2 3 4- s 2 4 + s 23 4- s 



24 



H, S^ 3 \ s 2 2 , Sg 4 , s 2 3 4- s 2 4 4- s 23 4- s 24 



35 



if, S^ 3 ), s 2 3 4- s 25 4- s| 5 , s 2 4 4- s 24 4- Sg 4 4- s 
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5 3 4 
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5 f 2 

4«-#3 
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o 
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yes 


• • 

O 5 


yes 


O 
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no 


si 
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• 

si 
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• 
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Table IA11 Spin graphs up to iV = 5 spins. 



